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The modulational instability of nonlinearly interacting spatially incoherent Stokes 
waves is analyzed. Starting from a pair of nonlinear Schrodinger equations, we derive 
a coupled set of wave-kinetic equations by using the Wigner transform technique. It 
is shown that the partial coherence of the interacting waves induces novel effects on 
the dynamics of crossing sea states. 

PACS numbers: 47.35.-i,47.15.Hg, 92.10.Hm, 02.60.Cb 



Extremely large amplitude 



lplitude 

fiflQB 



reak waves (also referred to as rogue or giant waves) in 



oceans are well-known |2|, y, IJ] . Freak waves, which are steep waves and which can appear 
quickly from a relatively calm responsible for the loss of many human lives and 

ships. Therefore, it is important to understand the nonlinear propagation q, y, q y, 3 of 

such waves in dispersive fluids that are far from a stationary state. Specifically, we recall 
that in non-stationary oceans, the dynamics of modulated Stokes waves is under certain 
circumstances governed by a nonlinear Schrodinger (NLS) equation [h3|, from which so 
called freak waves can occur. We note that the occurrence of freak waves may also be due 
to, e.g. a linear superposition of harmonic waves or wave-current interactions, and can thus 
be described by a different set of equations, such as the Zakharov equation The latter 
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depicts a modulational instability (also known as the Benjamin- Feir instability in 
fluid mechanics) of a constant amplitude carrier wave, and indicates that modula tonally 



unstable waves can form localized envelope disturbances JjJ. In Refs. [15( and 16( the 



effects of a random sea swell on the stability was investigated using the single sea state 
nonlinear deep water wave picture. 

Recently, Onorato et al. developed a two-dimensional weakly nonlinear model for 
two coherent Stokes waves in deep water with two different directions of propagation. For 



one-dimensional wave propagation, Ref. [17|] presented specific results for the modulational 
instability of two nonlinearly coupled coherent Stokes waves. It was found that freak waves 
may occur as a result of the modulational instability, and the growth rate for the crossing sea 
state modulational instability was found to be larger than that for the single sea state case. 
Since in reality the waves are incoherent it is however appropriate to examine the effects 
of random phases using perturbation theory. The presence of a spatial random phase of 
the background sea state gives rise to spectral broadening of the wave distribution function 
(as compared to the usual modulational perturbations around a monochromatic background 
wave distribution), and the modulational instability of nonlinearly coupled Stokes waves can 
be affected accordingly. This is the objective of the present Letter. 

In the following, we thus analyze the statistical properties of one-dimensional incoherent 
crossing see states. In particular, we investigate the modulational instability properties of 
these nonlinearly interacting waves. A generalized distribution function, which is valid for 
partially coherent waves, will be obtained. Using perturbation theory, it will be shown that 
partial coherence in terms of a random phase approximation gives rise to spectral broadening 
of the background crossing sea states. This broadening tends to stabilize the inherent system 
modulational instability, enabling the interaction between the waves over longer distances 
before any perturbation grows to appreciable levels forming freak waves. In fact, for a wide 
enough spectral distribution of waves the modulational instability can be suppressed. 



In Ref. [13] the equations 



and 



(1) 



3 



^-^ + ^| 2 + 2Cl^| 2 )5 = 0, (2) 

for one-dimensional comoving crossing sea states were presented. Here A and B represents 
the surface elevation of the crossing sea states, a is the group velocity dispersion parameter, 
£ is the nonlinear self-interaction parameter, and ( determines the strength of the nonlinear 
interaction between the crossing sea states. We note that the sign of a can always be 
chosen positive by a change of time coordinate. Thus, we take a > 0, while the nonlinear 
parameters may be positive as well as negative. It can also be noted that the one-dimensional 
nonlinear Schrodinger equation has been shown to be in excellent agreement with laboratory 
experiments [18 ]. 

There are traveling pure Stokes wave solutions A = A exp(— iuxt) and B = 
B exp(-iu B t) of Eqs. CQ) and ((21), where u A = £A 2 + 2(B% and u B = £5 2 + 2(A%. Letting 
A(t, x) = A [l + a(t, x)] exp(-iuAt) and B(t, x) = B [l+b(t, x)] exp(-iuBt), where a,b <^ 1, 
we linearize Eqs. (OQ) and (j2J), take the real and imaginary parts of the resultant equations, 



and Fourier analyze them against the 



The resulting dispersion relation is 



17| 



requency Q and wavenumber K of the modulations. 



fi 4 -2aifV^ 2 +^o+££o) fi2 +a 2 ^ 4 [("^ 2 + 2^)( a fr 2 + 2£5 2 ) - 16C 2 A 2 5 2 ] = 0, (3) 



i.e. 



1/2 



n = ± jaK 2 [aK 2 + Z(A 2 + B*)} ± aK^l^AlBl + ?(J% - B^j . (4) 

We now investigate the modulational instability of incoherent crossing sea states by ap- 
plying a Wigner transform to the wave amplitudes A and B. For a given function f(t, x), the 
corresponding Wigner function pAt, x, p) is defined as the Fourier transform of the two-point 



correlation function, i.e. 
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p f (t,x,p) 



1 

2^ 



dpe^(r{t,x + p/2)f(t,x-p/2)), 



(5) 



where the asterisk stands for the complex conjugate and the angular bracket denotes the 
ensemble average 26[. The function pf is a generalized distribution function for the waves 
of the sea states. The wave intensity If corresponding to the function / can then be written 

as 
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If=(\f\ 2 )= J dpp f (t,x,p). (6) 

Applying the time derivative to Eq. (jSJ), with / = A or B, and using the nonlinear 
Schrodinger equations (CQ) and (121), valid for slowly varying envelopes, respectively, we obtain 
the nonlinearly coupled wave-kinetic equations 



*' + ^-2( ft+ 2(y* iAilp^O, (7) 



dt dx " \ 2 dx dp 



and 



where the sin-operator is defined in terms of its Taylor expansion and the arrows give the 
direction of the differentiation. Equations (171) and (JSJ) model the nonlinear evolution of 
partially coherent crossing sea states. 

First order perturbations of flTJ) and ([S]) can be treated by letting pj(t,x,p) = Pjoip) + 
Pji(p) exp(iKx— iflt) and I j(t, x) = Ij +Ijiexp(iKx—iQt), where j = A, B, \pj\\ <C Pjo, and 
\Iji\ <C IjQ. Linearizing with respect to the perturbation variables, we obtain the dispersion 
relation from Eqs. (J7J) and (jHJ) 

(1 + £A A )(1 + £A B ) - 4( 2 A A A B = 0, (9) 

where 

3 J F n-2aKp v ; 

Here we have used 2z sin(iif/2S p )/)jo(p) = —Pjoip + ^/2) + PjoG — K/2). The dispersion 
relation ([9]) with ( flOl) is valid for partially coherent crossing sea states that interact nonlin- 
early. 

In the coherent case, the background waves are simple plane waves, so that the distri- 
bution functions are given by pj = Ij 5(p), where I A o = A% and I bo = Bq. Inserting the 
coherent background distribution function into (fTUl) . the dispersion relation (jUJ) reduces to 
Eq. @. 
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We now assume that the sea states suffer from random perturbations, such that they have 
partial phase coherence. Thus means that the phase (fj(x) of the background wave ampli- 
tudes A, B satisfies (exp[— i(fj(x + u/2) exp[i(pj(x — /x/2)]]) = exp(— Pjw\^\), corresponding 



to the Lorentz distribution function 23] 



where pjw represents the width of the jth distribution function. Thus, we see that the 
partial coherence in the waves phases give rises to a spectral broadening in terms of the 
distribution function. With this spectral background of sea states, the dispersion relation 
(l9l) takes the form 



2jaK 2 I m 
(Q + 2taKp BW ) 2 -a 2 K\ 

[(Q + 2iaKp AW ) 2 -a 2 K 4 ][(n + 2iaKp BW y -a 2 K 4 } ' 1 ' 
We next analyze the dispersion relation (fT2|) by letting, for simplicity, paw — Pbw — Pw- 
We then obtain 



2£aK 2 I A0 
CD, + 2iaKp AW ) 2 - a 2 K 4 



Q = -2iaKp w ± \aK 2 [aK 2 + £{I M + I B0 )] ± aK 2 ^/l6( 2 1 A0 I B0 + ^^ao ~ Ibo) 2 } , 

(13) 

which agrees with (jl]) when pw — > 0. We may then transform to dimensionless variables by 
Q -> Sl/y/a\£\, aK 2 aK 2 / y/a\£\, aKp w aKp w /^/a\^\, I j0 -> I j0 y/\^\/a, and ( -> 
(/£,. With these re-scalings, the group velocity dispersion a = 1, while the self-nonlinearity 
coefficient becomes £ — > ±1. In Figs. 1 and 2 we have displayed the dimensionless growth 
rate T = —iQ using Eq. ffT3|) . In Fig. 1 we have chosen I A0 = I B0 = 0.5, while in Fig. 2 we 
have I A0 = 51 B0 = 0.5. The lowering of the growth rate due to the spectral broadening of the 
background sea states pjo can clearly be seen. Thus, the effect of the partial coherence of the 
background sea states is to stabilize their dynamics, allowing for long distance propagation 
of envelope solitons without the growth of freak waves. It should be noted that in general the 
spectral broadening enables interaction between the ocean waves over longer distances before 
any initial perturbation grows to an appreciable size. However, for large enough spectral 
width p w the modulational instability growth rate may be completely suppressed. The 
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FIG. 1: The effect of partial coherence when a is normalized to unity and Iao = I bo = 0.5. In all 
the panels, the coherent (full) and incoherent (dashed) cases are compared. In (a) £ = £ = 1, in (b) 
£ = — £ = 1, in (c) £ = — £ = — 1, and in (d) £ = £ = — 1. The partial coherence stabilizes the wave 
modulation, i.e. the wave will be able to interact over longer distances before the perturbations 
grow to appreciable levels forming freak waves. In fact, for a wide enough spectral distribution, 
the modulational instability can be completely suppressed. 



results for pw = 0, i.e. for coherent background waves without spatial spectral broadening, 



are consistent with the results presented by Onorato, Osborne, and Serio [17| (see the full 
curves in Figs. 1 and 2), where it was found that crossing sea states have a larger growth 
rate for the instability as compared to the single sea state instability. 

Indeed, there are soliton solutions to the system ([I]) and Following Ref. [28j, we can 
find the dark-bright soliton pair 



A(t, x) = A tanh[X(t, x)] exp(inx — i5 1) 



(14) 
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FIG. 2: The effect of partial coherence when a is normalized to unity and Iao = 5/bo = 0.5. In all 
the panels, the coherent (full) and incoherent (dashed) cases are compared. In (a) ( = £ = 1, in 
(b) C, = — £ = 1, in (c) C = — £ = —1, and in (d) ( = £ = —1. As in the cases in Fig. 1, the partial 
coherence stabilizes the sea state modulation. 



and 

B{t,x) = B sech[X(t, x)} exp(inx — i5t), (15) 

where X = (x — Vt)/L, V = 2an, = B%, and 5 = an 2 + ^A^. Here A and L are treated 
as two free parameters determining the elevation and width of the solitons, respectively. 
Of course, the classical nonlinear Schrodinger equation has both soliton and multi-soliton 
solutions. Here we see that also the co-existence of dark and bright soliton water waves is 
possible. 

To summarize, we have presented an investigation of the modulational instability of two 
incoherent crossing sea states that are nonlinearly interacting in deep water. For this pur- 
pose, we have introduced the Wigner transformation of the coupled nonlinear Schrodinger 
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equations of Ref. [17| and obtained two coupled wave-kinetic equations (or von Neumann 
equations). The latter have been analyzed to obtain a nonlinear dispersion relation for back- 
ground sea states that have broadband spectra, i.e. finite spectral width. It is found that 
the growth rate of the modulational instability is suppressed. Hence, random phased non- 
linearly interacting waves could propagate over long distances without being much affected 
by the modulational instability, and for a wide enough spectral distribution the formation of 
freak waves is thus completely suppressed. However, it should be stressed that this complete 
suppression is a result of the above NLS model calculations. As an alternative, one could 
start from the original equations, and then consider the interaction of wave packets which 
are eigenfunctions of the corresponding linearised problem with a background of Stokes wave 
trains. In that case, a random mixture of the eigenfunctions will postpone, although not 
completely suppress, the freak wave appearance. 
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